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Using an approach inspired from Spin Glasses, we show that the multimode disordered Dicke 
model is equivalent to a quantum Hopfield network. We propose variational ground states for the 
system at zero temperature, which we conjecture to be exact in the thermodynamic limit. These 
ground states contain the information on the disordered qubit-photon couplings. These results lead 
to two intriguing physical implications. First, once the qubit-photon couplings can be engineered, it 
should be possible to build scalable pattern-storing systems whose dynamics is governed by quantum 
laws. Second, we argue with an example how such Dicke quantum simulators might be used as a 
solver of “hard” combinatorial optimization problems. 

PACS numbers: 


The connection of experimentally realizable quantum 
systems with computation contains promising perspec¬ 
tives from both the fundamental and the technological 
viewpoint [HE]. For example, quantum computational 
capabilities can be implemented by “quantum gates” [3] 
and by the so-called “adiabatic quantum optimization” 
technique BHB|. Today’s experimental technology of 
highly controllable quantum simulators, recently used for 
testing theoretical predictions in a wide range of areas of 
physics [Ml, offers new opportunities for exploring com¬ 
puting power for quantum systems. 

In the case of light-matter interaction at the quantum 
level, the reference benchmark is the Dicke model ! TU] . 
Studies of its equilibrium properties have predicted a su¬ 
perradiant transition to occur in the strong coupling and 
low temperature regime EM- The superradiant phase 
is characterized by a macroscopic number of atoms in the 
excited state whose collective behaviour produces an en¬ 
hancement of spontaneous emission (proportional to the 
number of cooperating atoms in the sample). Crucially, 
this phenomenology is in direct link with experimen¬ 
tally feasible quantum simulators. Recently, Nagy and 
coworkers [14] argued that the Dicke model effectively 
describes the self-organization phase transition of a Bose- 
Einstein condensate (BEC) in an optical cavity [1511 16| . 
Additionally, Dimer and colleagues m proposed a Cav¬ 
ity QED realization of the Dicke model based on cavity- 
mediated Raman transitions, closer in spirit to the orig¬ 
inal Dicke’s idea. Evidence of superradiance in this sys¬ 
tem is reported in [IS]. An implementation of generalized 
Dicke models in hybrid quantum systems has also been 
put forward m- More generally, Dicke-like Hamiltoni¬ 
ans describe a variety of physical systems, ranging from 
Circuit QED [20U24] to Cavity QED with Dirac fermions 
in graphene [25ff27l . Additionally, disorder and frustra¬ 
tion of the atom-photon couplings have an important role 



FIG. 1: In the Dicke model, photons (yellow lines) mediate 
a long range interaction between qubits (green circles). The 
drawing sketches schematically a six qubits system within its 
fully-connected graph and its internal level structure. In the 
standard single-mode Dicke model the exchange coupling is 
fixed at the same value for every pair of qubits. In systems 
where both many modes and disorder are present, the ex¬ 
change couplings are qubit-dependent and take the form given 
by Eq. Q. 

in the study of BEC in multimode cavities (28] [29]. Re¬ 
cent works (30, OH] discussed a multimodal-Cavity QED 
simulator with disordered interactions. The authors ar¬ 
gue that these systems could be employed to explore spin- 
glass properties at the quantum level [301 - i32] . However, 
the possible quantum computation applications of this 
new class of quantum simulators remain relatively unex¬ 
plored. 

In this Letter, we consider a multimode disordered 
Dicke model with finite number of modes. We calcu¬ 
late exactly (in the thermodynamic limit) the free en¬ 
ergy of the system at temperature T = 1//3 and we 
find a superradiant phase transition characterized by the 
same free-energy landscape of the Hopfield model [33] in 
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the so-called “symmetry broken” phase, with the typi¬ 
cal strong-coupling threshold of the Dicke model. From 
the theoretical standpoint, our results generalize to the 
case of quenched disordered couplings the remarkable 
analysis performed by Lieb et al. HM1- The choice 
of frozen couplings is compatible with the characteris¬ 
tic time scales involved in light-matter interactions. The 
calculation of the partition function leads us to suggest 
variational ground states for the model, which we con¬ 
jecture to be exact in the thermodynamic limit. 

The physical consequences of this analysis are fasci¬ 
nating: once the multimode strong-coupling regime is 
reached and qubit-photon couplings are engineered, it 
should be possible to build a pattern-storing system 
whose underlying dynamics is fully governed by quan¬ 
tum laws. Moreover, Dicke quantum simulators here ana¬ 
lyzed may be suitable to implement specific optimization 
problems, in the spirit of adiabatic quantum computa¬ 
tion mg. We point out a non-polynomial optimization 
problem mum, number partitioning, which could be 
implemented in a single mode cavity QED setup with 
controllable disorder. Computing applications based on 
cavity mediated interactions might owns the advantage 
to be a viable way to generate entagled many-body states 
with remarkable scalability properties, as recently shown 
in Ref. [35]. 

Hopfield’s main idea [33] is that the retrieval of stored 
information, such as memory patterns, may emerge as a 
collective dynamical property of microscopic constituents 
(“neurons”) whose interconnections (“synapses”) are re¬ 
inforced or weakened through a training phase (e.g. Heb- 
bian learning [361 EZ])- This is achieved in his model 
through a fictitious neuronal-dynamics whose effect is to 
minimize the Lyapunov cost function: 

1 N 

E=--J2 T vS i S j , (1) 

i,j— 1 

where N is the number of neurons, Si = 1 if the i-th neu¬ 
ron is active, and —1 otherwise, and the p stored patterns 
= ±1 (k = 1, • • • ,p) determine the interconnections 
Tij through the relation: T t j = 1 
The analysis in Ref. [33] shows that the long-time dy¬ 
namics always converges to one of the p stored patterns, 
i.e. these configurations are the global minima of the cost 
function 0. The interpretation of this result is that a 
suitable choice of the interconnections allows to store a 
given number of memory patterns into the neural net¬ 
work. Data retrieval is achieved through an algorithm 
that minimizes the energy function ([!]). A phase transi¬ 
tion to a “complex” phase marks the intrinsic limitation 
on the number of patterns p that can be stored. If p ex¬ 
ceeds the critical threshold p ~ 0.147V many failures in 
the process of retrieval occur [3SJ GUj. 

In this manuscript we consider the following multi¬ 


mode Dicke Hamiltonian: 

M N N,M 

h = '52 + A a i + 5Z 9ik (a k + 4) erf , (2) 

k — 1 2—1 i,k= 1 


effectively modelling quantum light-matter interaction 
of N two-level systems with detuning A and M elec¬ 
tromagnetic modes supposed to be quasi-degenerate at 
the common frequency u> and with couplings that we 
parametrize for future convenience as gi k = flgik/VN, 
where D is the Rabi frequency and the dimensionless gi k ’s 
are both atom and mode-dependent. In Cavity QED re¬ 
alizations, w represents the detuning between the cavity 
frequency and the pumping frequency and could be both 
positive or negative. A possible choice of the couplings 
is gu = cos (kiXi), being ki the wave vector of the photon 
and Xi the position of i-th atom [31:. 

We are interested in the thermodynamic properties of 
this system in the limit M <C N, and thus in evaluat¬ 
ing the partition function Z = Tve~^ H . This evaluation 
can be performed rigorously in the thermodynamic limit 
(N —> oo) using the techniques introduced in Refs. [TTL 
H3]. We first consider the fully-commuting limit A = 0. 
In this case the evaluation of the partition function is 
straightforward (see Supplementary material) and we ob¬ 
tain Z = ZfbZh, where Zfb is a free boson partition 
function and Zh is a classical Ising model with local 
quenched exchange interactions of the form: 


02 iV1 

j 9 ik 9jk 

Jij ~ 2 -^ 7 “ 

k =1 


(3) 


The physical interpretation of this result is that photons 
mediate long range interactions among the atoms, re¬ 
sulting in an atomic effective Hamiltonian described by 
a fully-connected Ising model (see Fig. [l]). The role of 
the couplings gik can be understood from Eq. (|T|) in the 
context of the Hopfield network. They are the memory 
pattern stored in the system. By computing exactly the 
free energy of the model, we will show that this inter¬ 
pretation stays unaltered in the more complicated case 
A ^ 0. 

We now proceed to the evaluation of the quantum par¬ 
tition function. We use the method of Wang and cowork¬ 
ers 1131 3D] (proved to be exact in the thermodynamic 
limit for M/N —> 0 [12]). We introduce a set of coherent 
states \a k ) with a k = x k + iy k , one for each electromag¬ 
netic mode k , and we expand the partition function on 
this overcomplete basis: 


„ M ,2 

z = I] — T D4<M|e-^|M>, (4) 

J k =1 ^ 

where Tr^ is the atomic trace only. The only techni¬ 
cal complication is the calculation of the matrix element 
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in Q. This turns out to be equal, apart from non- 
extensive contributions, to the exponential of the oper¬ 
ator in Eq ([2]) with the replacements a*,, a' k —► afe, a% 
mm■ At this stage the trace over the atomic degrees 
of freedom can be easily performed. The integral over 
the imaginary parts of a^’s give an overall unimportant 
constant. Finally, defining the M-dimensional vectors 
x = (xi,X 2 , ■ ■ ■ ,x M ) and g, = {gn, g i2 , ■ ■ ■ ,gm), and 
with the change of variables m = x/y/N, the partition 
function assumes a suitable form for performing a saddle- 
point integration, i.e. Z = f d M me~ N ^ m \ Here / is 
the free energy 


/(m) = /3m ■ m 


1 

N 


N 

^ log G(m, gj), 
1=1 


(5) 


with: G(m,gj) = 2 cosh 


/3 (A 2 + H 2 (gi • m) 2 ) 2 


The order parameter m describes the superradiant 
phase transition. Physically, it gives the mean number 
of photons in every mode m- Its value is determined 
by minimizing the free energy in Eq. I©. Solutions of 
this optimization problem are, in principle, gj-dependent, 
but in the thermodynamic limit both the free energy and 
the saddle-point equation are self-averaging [38]. Thus 
we conclude that the free energy and the saddle point 
equations are given by 


/(m) = /3m • m (log G(m, g)) g , 

m = T ( ^(gO S tanh 

1 

with: /r(g) = (A 2 + fl 2 (g • m) 2 ) 2 and (■ • ■) represent¬ 
ing the average over the disorder distribution. Eq. © re¬ 
duces to the mean-field equations for the Hopfield model 
for A —> 0 [38]. Thus, A may be intended as a quan¬ 
tum annealer parameter. To fully specify the model, the 
probability distribution for the couplings is needed. In 
the following we will assume 


p (g) = fl Q* (9k - 1) + (. 9k + 1)) , (7) 


but we have verified that the results are qualitatively ro¬ 
bust as long as the disorder is not too peaked around 
zero in accordance with the classical results of Ref. [35] 
To locate the critical point it suffices to expand in Tay¬ 
lor series Eqs. ©• As in the conventional Dicke model, 
a temperature-independent threshold fl 2 = 2A emerges. 
For fl < fl c , the phase transition is inhibited at all tem¬ 
peratures. Whenever the magnitude of the coupling ex¬ 
ceeds this threshold value, the critical temperature is lo¬ 
cated at T c = A/ arctanh (2A/H 2 ). 

Above the critical temperature T c the only solution to 
© is a paramagnetic state, with rrik = 0 for all k. Below 
T c , different solutions appear. We now set out to classify 



FIG. 2: Order parameter rrd 1 ' in Eq. © as a function of A in 
the ultrastrong coupling regime at T = 0. On the theoretical 
curve (dashed line) are superposed the results (green dots) 
of a naive numerical optimization algorithm minimizing the 
ground state energy ansatz in Eq. (111. The parameters of 
the simulation are fl = 4, N = 100, M = 10. Green dots are 
the result of a single disorder realization. This implies the 
correctness of our self-averaging hypothesis, which is a good 
approximation already for N = 100. The classical Hopfield 
model is recovered in the limit A —> 0. 


these solutions and their stability under temperature de¬ 
crease. For this analysis, we both considered the Hessian 
matrix d 2 f /drrikdmi (see Supplementary materials for 
its explicit expression) and numerical optimization (Fig¬ 
ure [2]) . The key point, as mentioned above is that in 
this “symmetry-broken” phase the system takes 2 M de¬ 
generate ground states (as well as many metastable states 
energetically well separated from the ground states). In 
other words, also in this fully quantum limit the free- 
energy landscape still closely resembles that of the Hop- 
field model [55] . 

The ground state solutions have the explicit form: 

m fe =m«(0,(),••• ,0, 1,0, - - - ,0 ). (8) 

v ^ ^ v 

k— 1 times M— k-\-l times 

Equation (|6j) for the order parameter m ^ reduces 
to: 2/u(m( 1 7) = H 2 tanh (/3/Lt(m^^)), where g,(m^) = 
•\/A 2 + Q 2 . In the zero temperature limit the or¬ 
der parameter can be evaluated exactly: 


J n2 

A 2 

V 4 

fl 2 


At zero temperature the most interesting state is the 
ground state (GS) of the Hamiltonian ©. Inspired by 
the calculation above we propose the variational ansatz 
for the GS: 


I GS) = |ai, a 2 , ■ ■ ■ , a M ) Ispin/aq, • • • , a M )) , (10) 

where |ai,a 2 j" -- ,«m) is the product of M coherent 
states and the spin part is factorized. The mean value of 
the energy in this GS is given by: 

Egs ( m ) = m ■ m — ^ \/A 2 + H 2 (g ■ m) 2 ^ . (11) 
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This expression exactly equals the free-energy computed 
previously in the limit f3 —> oo, which leads us to con¬ 
jecture that our factorized variational ansatz is exact. 
The quantum phase transition is located at the critical 
coupling = fl c = v / 2A, at which the paramagnetic so¬ 
lution becomes unstable. In the symmetry-broken phase 
we have 2 M degenerate ground states of the form 

I GS) k = | 0A ■ , 0 , ±m (1) , 0, • • • , 0 ) |spin(±m (1) )) fc , 

k—1 times M—k times 

( 12 ) 

with k = 1, • • • M. The spin wave function is also factor¬ 
ized |spin(±m (1) )) fc = n.; I Si) k , with 


controllable disorder as follows. Number partitioning can 
be formulated as an optimization problem [44]: given 
a set A = {ai, 02 , ■ • •, ajv} of positive numbers, find a 
partition, i.e. a subset A' C A, such that the residue: 
E = I E a^A' a i - E a&A' a J I is minim i z ed. A partition 
can be defined by numbers Sj = ±1: Sj = 1 if aj £ A ', 
Sj = —1 otherwise. The cost function can be replaced 
by a classical spin hamiltonian, whose ground state is 
equivalent to the minimum partition: 

N 

H = aiCLjSiSj . (14) 

i,j =1 


S i)k ~ 


1 ( -A + V / A ^+W k 


U 


Pik 


°i) T \di) 


(13) 


where A/” a normalization, /3 \k = gikrrP 1 ^ and |e*), \gp are 
erf’s eigenstates. It is worth noting that, as expected, the 
ground state energy is a self-averaging quantity, whereas 
the ground states are not, being disorder-dependent also 
in the thermodynamic limit. 

The above calculation shows that in the superradiant 
phase the ground state of the system is a quantum su¬ 
perposition of the 2 M degenerate eigenvectors given by 
Eqs. ppl ). Their explicit form suggests that at fixed 
disorder and mode number the information about the dis¬ 
ordered couplings belonging to the fc-th mode is printed 
on the atomic wave function. Moreover, the photonic 
parts of the wave functions are all orthogonal for k\ ^ k 2 
in the thermodynamic limit. This implies that in princi¬ 
ple a suitable measure on the photons-subsystem causes 
the collapse over one of the 2 M ground states and gives 
thus the possibility to retrieve information (“patterns”) 
stored in the atomic wave function. As mentioned above, 
a single-mode Dicke model has been recently realized 
with cavity-mediated Raman transitions in cavity QED 
with ultracold atoms jTS]. A Multimode cavity QED 
setup supporting disordered couplings has been proposed 
in refs. PM HQ, and preliminary evidence of superradi¬ 
ance in this system is found in j32j. A setup operating 
in multimode regime has been recently suggested also in 
Circuit QED [43] • We are not aware of a setup (in Cavity 
or Circuit QED) that implements both multimode strong 
coupling regime and controllable disorder, a necessary 
condition for the quantum pattern-retrieval system that 
we conjecture here. 

We surmise that Multimode Dicke quantum setups 
with controllable disorder could be used beyond stor¬ 
age, to simulate specific optimization problems. In¬ 
deed, finding the ground state of classical spin models 
with disordered interactions is equivalent, in most cases, 
to finding solutions of computationally expensive non¬ 
polynomial (NP) problems [51] , For example, the sim¬ 
plest NP-hard problem, number partitioning, could be 
implemented in a single-mode cavity QED setup with 


In a single mode cavity QED network couplings can 
be chosen as gi = cos(kxi) pH]. By the definition of 
a = max _4 aj and a,j = aj /a, it is possible to engineer 
the gP s in such a way to implement a given instance of 
the problem provided that the cavity is in the “blue” 
detuned regime to ensure the appropriate sign for the 
couplings, see Eq. § With a suitable adiabatic anneal¬ 
ing of the atomic detuning A, the system should collapse 
on qubit configurations that are good solutions of the 
corresponding optimization problem. 

In conclusion, this Letter provides the first rigorous 
analysis of the multimode disordered Dicke model, valid 
beyond the weak-coupling regime and exact in the ther¬ 
modynamic limit. The equivalence between multimodal 
disordered Dicke model and a quantum Hopfield net¬ 
work [45], together with the proposal of a cavity QED 
setup implementing a non polynomial optimization prob¬ 
lem, demonstrates the possibility of quantum computa¬ 
tional abilities of this new class of quantum simulators. 
Our proposal is conceptually complementary to a stan¬ 
dard quantum computation perspective mm- Indeed, 
the information can be “written” on the qubits through a 
quantum annealing on the detuning A, similarly to what 
happens for adiabatic quantum computation BHS|. 
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supporting material 

Pietro Rotondo, Marco Cosentino Lagomarsino, and Giovanni Viola 

In this material, we give more details on the derivations of the results presented in the main text. 


DERIVATION OF EQUATION (3) 


We begin with the derivation of Eq. (3), of the main text. We consider the partition function Z = Tve~^ H with H 
given in Eq. (2) for A = 0. In this fully commuting limit we can evaluate the partition function straightforwardly. 
We introduce new set of bosonic operators: 


bl = at 


n 0 

/ y 9ik&i ) bk = dk H 7 = / y 9ik&i ■ 


1=1 


1=1 


with [bk 1 , &£,] = <5fc,fe'. By means of those, H is written as the sum of two commuting operators: 


(S.l) 


M „2 N M 

h =u>j 2 b i b k - g,i.g J k<'T' l (s.2) 

k =1 i,j =1 k=l 

As a byproduct we obtain the factorization of the full partition function ISlj :Z = Zg Zg 1 where Zg is an overall free 
boson partition function that we can safely ignore in the thermodynamic limit. On the other hand 


Zg = Tr^exp /3 y ^ JjjOjOj 


(S.3) 




is an Ising contribution with both spin and mode dependent couplings of the form given in Eq. (3) of the main 
text. In Eq.(S.3) Tr CT indicates the trace on the spins only. 


DERIVATION OF EQUATION (5) 


In this section we report the derivation of Eq. (5) of the main text, which essentially follows the derivation of Wang 
and Hioe [S2] proved to be rigorous by Hepp and Lieb [S3| . The authors of Ref. [52] have shown explicitly, in the 
termodynamic limit, that the convenient way to calculate the trace on the Hilbert space of bosons, in the partition 
function, is to evaluate it on the set of the coherent states |{a}). The photonic matrix element in the partition 
function of Eq. (4) equals in the thermodynamic limit (a*, = Xk + iyk)' 


M 


N 


n 


N M 


(M|e pH IM) ^ exp ( -pu^2(xl + y 2 k ) - ~ 


\ k=l 

The atomic trace thus factorizes and it can be calculated: 


Qik^k^i 


(S.4) 


i =1 k= 1 


Z = f [] e ~^^l + vl) J^ C osh L* A + ^(f 

J k=l ^ 2 = 1 \ ) \fc = l 


9ik%k 


, M 

n 

k=1 


dm k Nf(mj3) 


(S.5) 


In the last term we introduced the vectorial notation defined in the main text. The final expression for the free energy 
is: 


1 N 

- ^yE iogc ° sh 



/(m, ff) = /3m • m 


(S.6) 
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By minimizing the free energy above and using the self averaging property of Eq. 
field equations: 


(S.6), we obtain the exact mean 


Q 2 


rn = 




To locate the critical point it suffices to expand in Taylor series Eqs. (S.6),(S.7): 

n 2 


/(m) - /(0) =/3[l- — tanh(/3A) ) m • m + O ( m \) 


Q 2 

m k = — tanh(/3A) m k + 0(r 


(S.7) 


(S.8) 


As in the conventional Dicke model, a temperature-independent threshold fl 2 = 2A emerges. For fl < fi c , the phase 
transition is inhibited at all temperatures. Whenever the magnitude of the cou plin g exceeds this threshold value, the 
critical temperature is located at T c = A/arctanh (2A/fi 2 ). Solutions to Eq. (S.7) can be classified according to the 
number of non-zero components n of the order parameter m IS4| : 


m„ = m 


(n) 


( 1 , 1 ,--- , 1 , 0 , 0 ,--- , 0 ), 


(S.9) 


n times M—n times 


where all permutations are also possible. In particular, solutions for n = 1 are the ones with the lowest free energy. 
There are 2 M of such solutions, corresponding to the Z 2 x Sm symmetry breaking of our multimode Dicke model. 

For completeness we report the explicit expression for the Hessian matrix of the free energy, omitted in the main 
text for space imitations: 


d 2 f 

dm k dmi 


2Su - n 2 tanh(/3/x(g))^ + H 4 


(g ■ m) 2 g k gi 

Me) 2 


tanh(/3/x(g)) 

Mg) 


/3sech 2 (MMg)) 


(S.10) 


[51] P. Rotondo, E. Tesio and S. Caracciolo, Phys. Rev. B 91 014415 . 

[52] K. Wang and F. T. Hioe, Phys. Rev. A 7 831. 

[53] K. Hepp and E. Lieb, Phys. Rev. A 8, 2517. 

[54] D. J. Amit, H. Gutfreund and H. Sompolinsky, Phys. Rev. A 32 1007. 
















